Observation of the spontaneous vortex phase in the weakly ferromagnetic 
superconductor ErNi2B2C: A penetration depth study 
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The coexistence of weak ferromagnetism and superconductivity in ErNi2B2C suggests the possi- 
bility of a spontaneous vortex phase (SVP) in which vortices appear in the absence of an external 
field. We report evidence for the long-sought SVP from the in-plane magnetic penetration depth 
AA(T) of high-quality single crystals of ErNi2B2C. In addition to expected features at the Neel 
temperature Tn = 6.0 K and weak ferromagnetic onset at Twfm = 2.3 K, AA(r) rises to a maxi- 
mum at Tm = 0.45 K before dropping sharply down to ~0.1 K. We assign the 0.45 K-maximum to 
the proliferation and freezing of spontaneous vortices. A model proposed by Koshelev and Vinokur 
explains the increasing AA(T') as a consequence of increasing vortex density, and its subsequent 
decrease below Tm as defect pinning suppresses vortex hopping. 



It is now clear that the borocarbide superconductor 
ErNi2B2C develops weak ferromagnetism (WFM) below 
"^WFM = 2.3 K while remaining a singlet supercon- 
ductor P, 0. The question naturally arises: how do 
these two seemingly incompatible orders — ferromag- 
netism and superconductivity — coexist microscopically? 
Clearly superconductivity will be suppressed if the in- 
ternal field Bin generated by the ferromagnetic moment 
exceeds He for a Type-I, or Hc2 for a Type-II, super- 
conductor (SC). For a Type-II SC, however, vortices are 
predicted to appear spontaneously if Bin lies in the range 
Hci < Bin ~ 47rM < ffc2 HllllH- In this sponta- 
neous vortex phase (SVP), the vortex screening currents 
shield superconducting regions from the intrinsic mag- 
netization. The vortices, however, may be qualitatively 
different from those generated by externally applied fields 
. In this Letter we report unusual features in the pen- 
etration depth data of a high quality single-crystal of 
ErNi2B2C that give strong evidence for the existence of 
the SVP. 

There have been previous SVP reports that we con- 
sider inconclusive. Ng and Varma |8|, for example, in- 
terpreted small angle neutron scattering (SANS) data on 
ErNi2B2C as a prelude to the SVP. In that experiment, 
Yaron et al. reported that the vortex-line lattice be- 
gins to tilt away from the c-axis (along which the mag- 
netic is field applied) towards the a-b plane below Twfm- 
However, the tilt can merely be a result of the vector sum 
of the applied field and the internal field produced by the 
ferromagnetic domains in the basal plane. Additional ev- 



idence was provided by SANS data with the applied 
magnetic field in the basal plane. A large field was ap- 
plied to align ferromagnetic domains. When the field was 
removed, the flux line lattice was found to persist below 
Twfm but disappear above it. However, owing to the 
low Tc (~ 8.5 K) and increased pinning below Twfm 
trapped flux cannot be ruled out. 

Among the magnetic members of the rare-earth (RE) 
nickel borocarbide family, RENi2B2C (RE — Ho, Er, 
Dy, etc.), ErNi2B2C, is a particularly good candidate 
for study. Superconductivity arises at Tc ~ 11 K and 
persists when antiferromagnetic (AF) order sets in at Tn 

6 K ■ In the AF state the Er spins are directed along 
the &-axis, forming a transversely polarized, incommen- 
surate sinusoidal spin-density-wave (SDW) state, with 
modulation vector modulation vector 6 = 0.553a* (a* 
The appearance of higher-order reflections 
at lower temperatures Q signals the development of a 
square-wave modulation, with regular spin slips spaced 
by 20a. Below 2.3 K WFM appears with B,n =0.1 T, 
approximately one Er magnetic moment per twenty unit 
cells, clearly correlated with spin slips. 

The relative stability of various phases of a ferro- 
magnetic superconductor was explored by Greenside 
et al. A spiral phase is not possible in the presence 
of strong uniaxial anisotropy and the spontaneous vor- 
tex phase is more stable than a linearly polarized state 
for small values of C = [Ffm/Fs], the ratio of fer- 
romagnetic to superconducting free-energy densities at 
T = 0. For C = 100 and the ratio A/7 = 10, where 
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7 = [?,kBT^S/{2aM'^{S + l)]^/^ is related to the ex- 
change stiffness, Greenside et al. find the SVP to be 
the most stable low-temperature phase; indeed, they sug- 
gest that the effect is most likely to be found in a dilute 
ferromagnetic superconductor, and that smaller values 
of C favor SVP. In the case of ErNi2B2C, where only 
5% of the Er atoms contribute to ferromagnetism, we 
have Fs — ~HI/%tt k -1.5 x 10^ erg/cm^, where Hc~ 
1900 G from Ref. The ferromagnetic energy density 
is Ffm = -3NkBTcS/{2{S +1)) w -4.3 x 10^ erg/cm^, 
where N— 1.5x10^^ cni~^ is the density of the (magnetic) 
Er atoms, and 5* = 3/2 is the Er spin. This then gives 
C = 30, strongly favoring the SVP. The spin-stiffness 
length is 7 = 100 A at low temperatures where M « 88 
G, so that A/7 « 7, close to the value assumed in Ref. |3 
As ErNi2B2C is strongly Type II (A/^ w 5), we conclude 
that the SVP phase is the preferred state for coexisting 
ferromagnetism and superconductivity. 

We have measured the temperature dependence of the 
in-plane magnetic penetration depth AA(r) = A(T) — 
A(Tbase), in single crystals of ErNi2B2C down to Ttase = 
0.12 K using a tunnel-diode based, self-inductive tech- 
nique at 21 MHz [3| with a noise level of 2 parts in 
10^ and low drift. The magnitude of the ac field was 
estimated to be less than 40 mOc. The cryostat was sur- 
rounded by a bilayer Mumetal shield that reduced the dc 
field to less than 1 mOe. The very small values of the ac 
and dc field in our system ensure that our measurement is 
essentially a zero-field one, thereby eliminating the pos- 
sibility of trapped flux. Details of sample growth and 
characterization are described in Ref. |l^ The samples 
were then annealed according to conditions described in 
Ref. ^3 The sample was mounted, using a small amount 
of GE varnish, on a single crystal sapphire rod. The other 
end of the rod was thermally connected to the mixing 
chamber of an Oxford Kelvinox 25 dilution refrigerator. 
The sample temperature is monitored using a calibrated 
Ru02 resistor at low temperatures {Thase ~ 1-8 K), and 
a calibrated Cernox thermometer at higher temperatures 
(1.3 K - 12 K). 

The deviation AA(r) = A(r) - A(0.12 K) is propor- 
tional to the change in resonant frequency A/(T) of the 
oscillator, with the proportionality factor G dependent 
on sample and coil geometries. We determine G for a 
pure Al single crystal by fitting the Al data to extreme 
nonlocal expressions and then adjust for relative sample 
dimensions (IQ] . Testing this approach on a single crystal 
of Pb, we found good agreement with conventional BCS 
expressions. The value of G obtained this way has an 
uncertainty of ±10% because our sample, with approxi- 
mate dimensions 1.2 x 0.9 x 0.4 mm"^, has a rectangular, 
rather than square, basal area [l^. 

Figure n shows the temperature-dependence of the in- 
plane penetration depth AA(T). We see the following 
features: (1) onset of superconductivity at T* = 11.3 K, 
(2) a slight shoulder at Tn = 6.0 K, (3) a broad peak at 
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FIG. 1: Temperature dependence of the penetration depth 
Alambda{T) from 0.12 K to 13 K. Inset shows the low- 
temperature region. The arrows point to the features at Tjv 
and TwFM- 



TwFM — 2.3 K, (4) another sharp peak at Tm = 0.45 K, 
and (5) an eventual downturn below T,„. The features at 
Tn and Tw f m have not been seen in previous microwave 
measurements of AAfT) on either thin-film 18] or single- 
crystal ErNi9B'jG Jl9l| but the former has been observed 
in SANS data |20j. We show in a separate publication 
(21] that the feature at Tn is only observed for relatively 
small non-magnetic scattering rates. The large value of 
Tc and the resolvability of the features at T/v and T-^^fm 
attest to the high purity of the samples. We show, for 
comparison, data for a sample grown by floating-zone 
methods in the inset to Fig. [3 No clear signal is seen at 
T/v, although there may be some sign of the Neel transi- 
tion near 5 K. In place of the up-turn in the penetration 
depth, the signal levels off near Tw fm before decreasing 
below 1 K. This suggests that spontaneous vortices at 
the surface of this sample are strongly pinned. 

Figure 121 shows the data below Twfm- The strong up- 
turn is a significant deviation from the normal monotonic 
decrease of the penetration depth with decreasing tem- 
perature. Because we expect the Meissner effect to vanish 
(A(r) to diverge) in the SVP in the absence of pinning |^, 
it is natural to analyze the low-temperature data in the 
context of weakly pinned vortices in the low-frequency 
limit. We use a two- level tunneling model proposed by 
Koshelev and Vinokur f KV) |22l| . This approach has 
been revisited by Korshunov |23| | and applied to ultrathin 
cuprate films by Calame et al. j2^. At relatively high 
frequencies, small oscillations of the pinned lattice near 
equilibrium (Campbell regime) dominate absorption. At 
lower frequencies, jumps of lattice regions between dif- 
ferent metastable states (two-level systems) come into 
play and determine the absorption. Both regimes are 
sensitive to the pinning strength, which depends on the 
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FIG. 2: (O) Temperature dependence of the penetration 
depth AA(r) in the WFMphase. We've assumed X{Ttase) ~ 
A(0)=700 A from Ref. The sohd hne is the fit of the 

data to Eqn. for n—0.21, while the dotted hne is that 
of n=0.125. The values of other parameters are mentioned 
in the text. The inset shows AA for a floating-zone-grown 
sample exhibiting no signals at Tjv or Twfm 



value of the internal magnetic field Bin- At small field, 
vortices are pinned independent of one another. When 
Bin exceeds a characteristic value Bp, pinning becomes 
collective 0|, and the vortex lattice spfits into volumes 
that are pinned as a whole, with correlated regions hav- 
ing length Lco parallel to the field. We argue below that, 
in the WFM phase, Hci < Bin, Bp ^ Hc2 and that the 
frequency lies in the two-level regime. 

In the two-level model, the pinning state is character- 
ized by a large number of neighboring metastable config- 
urations. If we shift any volume V of the vortex lattice, 
then a finite probability exists that at some distance u 
there is another state with an energy within A of the 
starting configuration and separated from it by a barrier 
U . At finite temperatures such regions of the lattice jump 
among metastable states. Under the action of an applied 
ac field of frequency to, the ac screening current exerts 
an ac Lorentz force on the vortex lattice which induces 
jumps among nearby metastable states. The motion of 
the vortices, which carry magnetic fields, then increases 
the ac field penetration into the sample, resulting in an 
increase in the effective penetration depth Ae/ / . The sys- 
tem exhibits typical Debye-type behavior with the prop- 
erties determined by t(T) ~ tq exp(?7/r), the mean time 
between jumps. When two-level response dominates over 
Cambell behavior the penetration is estimated jl^ to be- 
have as 



Xlff{T) = Xl{T) + 



Bl{T)nti 
IGttT 



(l + (tJT(T))2)cosh^(A/r) 



two level systems. Here A||(T) is the London penetra- 
tion depth in the absence of vortices; Bin, the internal 
magnetic field; and riti, the concentration of two-level 
systems. In the low- field region {Bin < Bp), the vortex 
lines move independently, and their presence does not 
change the penetration depth considerably (Ag// ~ A||). 
However, in the collective pinning state {Bin > Bp), the 
jumping volume is not too small, and the characteristic 
distance at which the nearest metastable state exists is 
approximately the radius of the pinning force u « ^|| . 
As the temperature decreases, there is insufficient ther- 
mal energy to overcome the barrier U. No jumping takes 
place, the vortices are frozen, and hence there is no extra 
penetration. One therefore recovers the London penetra- 
tion depth A|| at the lowest temperatures. 

The solid line shows the fit of Eq. 2] to the data be- 
low Twfm- In this fit, we follow KV and replace (...) 
with values that characterize an effective number n^f^oi 
active two-level systems. The values of the following 
quantities will be justified later: B{T = 0) « 1100 G, 
w « ^11 « 150 A, V ^ Lcow^ = 5.4 X IQ-i^ cm^, and 
Tq = 2.2 X 10~^ s. The temperature-dependence of the 
internal magnetic field B{T) can be obtained by fitting 
magnetization values in Ref. |3 to the expression 



B{T) 



1 - 



T 



T, 



WFM 



(2) 



(1) 

where (...) denotes an average over the distribution of 



giving n=0.21. This value of n is between the 2D-Ising 
value of 0.125 and the 3D-Ising value of 0.31, which is 
reasonable because in ErNi2B2C the spins lie on sheets 
normal to the a axis and are confined to be along or 
anti-along the b axis, yet there is also 3D behavior in 
the superconductivity. Because U and A are strongly 
correlated, we make the reasonable assumption that all 
metastable states are equivalent (A = 0) and choose 
the energy barrier U — 0.49 K and pinning density 
n^J^ = 1.61x10" cm-3 that best fit the peak in Ae//(r). 
This value of the barrier makes cut « 1 near 1 K. Note 
that this value of U is close to the position of the peak 
at Tm — this is reasonable since below this tempera- 
ture, the vortices no longer have enough thermal energy 
to overcome the barrier to hop among metastable states; 
hence, one recovers the Meissner state with A decreas- 
ing. We expect A||(T) to exhibit a power-law tempera- 
ture dependence at low temperatures from the combina- 
tion of gap-minima observed in non-magnetic borocar- 
bides and the increased pair-breaking as Er spins dis- 
order. Consequently, we set A||(T) = \\{0){1 + bT^) 
with b = 0.036 the third adjustable parameter in 
the fit. For comparison, we show the (dotted-line) fit 
with n = 1/8 (2D-Ising model) for which U = 0.49 K, 
n^J^ = 1.57 X 10" cm-3, and b = 0.033 K^^. Both fits 
reproduced the qualitative features of the data, though 
ihe latter curve fits the data slightly better. 

To justify our application of the two-level model to 
our data, we evaluate various physical parameters in 
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Quantity Expression Value Notes 

Depairing current, j\, c'l>o/(12^/37^^^ll Ajj) 1.5 x 10* A/cm^ 

Viscous drag coefficient, r; ~ $Q/(27r^y p„c^) 5.6 x 10~^ erg s cm~'^ Bardeen-Stephens model 

Bean-model critical current, jc 4cMh/L 1.9 x 10* A/cm^ Ref. fTH . L — 1 mm, 

Mh from hysteresis loop 

Flux coherence length, Lco ^'^ ^ Ina;; x = 19.5 nm Ref. El, [23, Ax/Ay 1.3 

Jump time prefactor, to V^wc/i^ojc) 2.2 x 10"^ s Ref. [2l | 

Collective pinning field, Bp $o(lnx)''/V-^'co 20 Oe Ref. [13 

Campbell crossover, w^r{B,T) ^ T^o/{rjBVf) 27 MHz Ref. [2^, -8=500 Oe; T=2 K 

TABLE 1: Vortex parameters for the two-level hopping regime 



the model using standard expressions for the vortex 
state [23. We start virith the measured quantities: 
zero-temperature in-plane penetration depth A||(0) ~ 
700 A and coherence lengt h £ || « 150 A (and k — 
— 4.7) from Ref. |2fl|, ferromagnetic moment 
M « 0.62/XB/Er « 100 Oe well below Twfm from 
and the normal-state resistivity Pn{T*) = 5.8 ^Q, 
cm that we measured for our sample. The value of M 
corresponds to an internal field B ^ AttM — 1100 G, 
which is greater than Hd 500 G, putting us in the 
mixed state. 

In Table I, we give the expressions and values for the 



quantities that lead to the flux coherence length Lco (19.5 
nm), the collective pinning field Bp (20 Oe), the fre- 
quency LOcr above which Campbell response is expected 
(27 MHz), and the jump-time prefactor tq (2.2 ns). Since 
we operate at 21 MHz, this puts us in the two-level 
regime. Note that Bp is less than Hd, suggesting the 
the mixed state of ErNi2B2C is always in the collec- 
tive pinning regime. Based on these values, we esti- 
mate the maximum density of two-level volumes to be 
~ 2.4 X lO^'^cm"'^ at the lowest temperatures, indicat- 
ing that approximately 1% are active in our frequency 
window. 



KV [23 also found the Campbell penetration depth 
to be r^B^, which is a monotonically increasing function 
with decreasing temperature, i.e. there is no peak at low 
temperatures. This is in agreement with our not being in 
the Campbell regime. We also measured A(r) with the ac 
field along the basal plane, finding features qualitatively 
similar to the present data, including the strength and 
position of the features at Tpf, Twfm and T^,. 

In conclusion, penetration depth data of single-crystal 
ErNi2B2C down to ~0.1 K provide strong evidence 
for the existence of a spontaneous vortex phase below 
Tw F M ■ The high quality of our sample enables us to see 
features at Tj^ and that have not been observed in 

previous studies of the penetration depth [11 [13. Other 
samples, such as that shown in the inset to Fig. [3 show 
no clear signal at either Tjv or Twfm, nor the upturn 
in penetration depth that we attribute to weakly pinned 
spontaneous vortices. As pointed out by Radzhiovsky [3, 
the SVP lattice is much softer than a conventional lattice, 
and therefore especially sensitive to quenched disorder. It 
may well be that the spontaneous vortices may be glass- 
like rather than forming a lattice, and that ferromagnetic 
closure domains form at the surface. While these aspects 
may make spontaneous vortices difficult to detect by neu- 
tron scattering or surface magnetization, vortices will still 
have strong effects on the electrodynamics, as observed 



here. 
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